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Abstract 

Process  simulation  has  become  an  accepted  tool  for  the  performance,  design,  and  optimization  of 
thermal  processes.  Solving  the  mathematical  models  representing  solar  heating  process  units  and 
systems  is  one  of  the  most  tedious  and  repetitive  problems.  Nested  iterative  procedures  are  usually 
needed  to  solve  these  models.  To  tackle  these  problems,  several  researchers  have  developed  different 
methods,  techniques,  and  computer  programs  for  the  simulation  of  very  wide  verity  of  solar  heating 
process  units  and  systems. 

It  is  of  interest  in  this  work  to  characterize  and  classify  these  methods,  techniques,  and  programs  in 
order  to  better  understand  their  relations,  types,  structures,  and  procedures. 

The  simulation  problems  are  outlined;  the  simulation  programs  are  grouped  into  two  main  types; 
special  purpose,  and  general-purpose  programs.  Sequential  and  simultaneous  computational 
sequences  are  illustrated.  Simulator  structure,  program  evaluation,  and  numerical  techniques  are 
summarized. 

By  considering  the  unit  and/or  system  entropy  generation  as  well  as  the  energy  and  material 
balances  equations,  more  realistic  models  can  be  obtained.  Also,  rapid  development  of  computer 
hardware  and  software  will  suggest  new  techniques  and  programs  to  be  considered.  These  progress 
directions  are  noted. 
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1.  Introduction 

The  demand  for  efficient  solar  heating  processes  requires  suitable  computational  tools 
for  designing  and  evaluating  new  and  existing  processes.  Process  simulator  (also  called 
process  flow  sheeting)  has  become  an  accepted  tool  for  this  purpose.  Also,  it  is  currently 
used  in  both  industry  and  engineering  education.  Any  attempt  to  gain  information  about  a 
process  by  carrying  out  calculations  on  a  mathematical  model  of  that  process  is  defined  as 
a  simulation  [1,2]. 

The  mathematical  model  representing  material  and  energy  balance  of  a  solar  heating 
process  may  comprise  a  large  system  of  equations  characterized  by  complexity,  non¬ 
linearity,  and  implicitly.  The  nested  and  countercurrent  information  flow  nature  of  the 
problem  adds  to  the  difficulty  of  the  solution  of  such  process.  With  a  flexible  computer 
program,  a  large  number  of  flow  sheeting  problems  can  be  tackled.  These  problems  can  be 
generally  divided  into  three  classes:  (i)  performance  problems,  (ii)  design  problems,  and 
(iii)  optimization  problems. 

Many  ‘special’  computer  programs  designed  for  different  solar  heating  system 
calculations  have  been  reported  [3].  However,  when  only  ‘general’  programs  are 
considered,  i.e.  programs  not  restricted  to  a  certain  type  of  calculations,  to  a 
fixed  flow  sheet,  certain  physical  data,  only  a  few  programs  can  be  found  in  the 
literature  [4]. 
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Without  properly  classified  programs,  methods,  techniques,  and  problems,  one  cannot 
hope  to  make  headway  in  the  field  of  solar  heating  process  simulation.  In  fact,  no  such 
classification  is  available  in  the  open  literature  for  this  subject. 

The  main  purpose  of  this  article  is  to  present  an  overview  of  the  main  lines  of  the  design 
and  modeling  philosophy  utilized  in  developing  solar  heating  process  simulation 
programs.  A  classification  of  the  concepts,  programs,  and  techniques  of  a  large  number 
of  available  mathematical  models  for  these  processes  are  also  outlined.  Of  course,  because 
of  the  space  and  time  limitations  the  paper  cannot  go  in  details. 


2.  Simulation  problems 

Consider  a  typical  solar  heating  process,  as  shown  in  Fig.  1,  and  suppose  that  one  wants 
to  simulate  its  performance.  The  first  step  in  the  modeling  of  the  process  is  to  structure  it. 
This  leads  to  a  flow  sheet  of  coupled  blocks  that  can  be  modeled  independently,  Fig.  2. 
Each  of  lettered  blocks  represents  a  processing  unit.  Each  of  these  units  is  connected  to  at 
least  one  of  the  other  units  by  one  or  more  pipes  or  streams  as  indicated  by  the  numbered 
line.  Arrows  into  a  box  indicate  streams  that  must  be  known  in  order  to  simulate  that  unit. 
Arrows  out  of  the  box  indicate  streams,  which  are  computed  from  the  process  simulation 
mathematical  model 

Once  it  has  been  decided  which  streams  are  inputs  to  each  unit  and  which  are  outputs, 
some  connection  equations  between  the  two  must  be  determined  as;  output  variables  =  F 
(input  variables).  This  model  equation  of  the  individual  units  is,  of  course,  the  most 
important  aspect  of  process  simulation.  Generally,  the  input  variables  are  the  inlet  feed 
stream  variables  and  the  model  parameters.  The  output  variables  are  the  outlet  stream 
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Fig.  2.  Flow  diagram  for  solar  heating  process. 
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variables,  performance  and  sizing  results.  Also,  the  various  physical  parameters  such  as 
transfer  area,  transfer  coefficient,  etc.  are  contained  in  the  model  equations. 

Now,  if  the  process  flow  diagram  is  without  recycle  streams  and  the  feed  streams  are 
known,  it  is  then  possible  to  calculate  the  whole  process  by  working  through  the  units 
sequentially.  This  is  the  most  straightforward  form  of  the  simulation  problem.  However,  if 
the  process  does  have  one  or  more  recycle  streams  as  input,  one  must  make  initial 
estimations  of  all  these  recycle  streams.  This  leads  to  a  set  of  computed  values  for  each  of 
the  recycle  streams  that  are  then  compared  with  the  original  estimations.  Some  iterative 
procedure  is  then  used  to  continue  calculations  around  the  units  (in  the  recycle  loop)  until 
the  streams  match  each  other. 

Iterations  may  exist  at  several  levels,  for  instance,  stream  2  depends  not  only  on  the 
value  of  stream  1,  but  it  also  depends  on  the  value  of  stream  3.  Nevertheless,  if  we  tear 
(break)  streams  3  and  5  then  use  an  iterative  procedure  until  the  values  of  tom  variables  are 
equal  to  the  values  obtained  by  the  calculation,  then  the  value  of  stream  6  (product)  is 
obtained.  Another  level  of  iteration  may  arise  because  of  the  complexity  of  the 
mathematical  models  used  to  describe  the  process  units.  For  example,  if  the  heat  storage 
tank  (unit  C)  is  modeled  by  a  non-linear  mathematical  model.  This  model  will  be  solved 
iteratively  to  calculate  streams  5  and  6  form  streams  4  and  7.  Finally,  the  third  level  of 
iteration  may  arise  if  one  wants  to  calculate  the  feed  (stream  1)  in  order  to  obtain  a  product 
of  given  specification  (stream  6).  This  requires  iteration  on  the  values  of  streams  2-5,  and 
7  until  stream  6  has  the  specified  values.  So,  it  is  quite  clear  why  simulation  is  a  real 
problem.  Several  types  of  problems  can  be  solved  by  a  simulation  program  (flow  sheeting 
program)  [5]. 

In  the  performance  problem ,  the  variables  associated  with  the  feed  streams  to  a  process 
unit  and  all  design  parameters  (such  as  solar  collector  area,  heat  exchanger  area,  etc.)  are 
assumed  to  be  known.  The  variables  associated  the  internal  and  output  streams  are  the 
unknowns. 

In  the  design  problem ,  some  design  parameters  and/or  feed  variables  are  left 
unspecified  and  become  unknown.  A  corresponding  number  of  additional  equations 
(equality  constraints)  relating  some  of  these  variables  are  added,  such  that  the  total  number 
of  unknowns  equates  the  number  of  equations. 

The  optimization  problem  differs  from  the  design  problem  in  that  the  number  of 
equality  constraints  is  smaller  than  that  of  the  variables  left  unspecified.  The  unspecified 
variables  are  now  calculated  so  as  to  minimize  an  objective  function,  normally  of 
economic  nature.  In  this  case,  the  solution  can  be  considered  also  by  using  inequality 
constraints. 


3.  Simulation  programs 

Many  computer  programs  for  different  solar  heating  processes  have  been  reported  in 
the  literature  [3,4].  These  programs  can  be  grouped  into  two  main  types:  special  purpose 
programs,  and  general-purpose  programs. 
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3.1.  Special  purpose  programs  (one-off  programs) 

For  a  particular  solar  heating  process  (or  unit)  with  a  fixed  configuration,  it  is  possible 
to  write  a  mathematical  model  describing  that  process  is  in  the  form  of  a  ‘one-off’ 
computer  program.  The  structure  of  these  programs  is  rigid,  simple,  and  straightforward. 
All  that  the  user  has  to  supply,  is  the  data  and  the  executive  handles  the  program  in  the 
same  way,  irrespective  of  the  nature  of  the  process  simulated. 

The  disadvantage  of  such  programs  is  that  a  model  exists  for  only  one  process  and  any 
changes  made  to  that  process  might  require  extensive  re-programming.  However,  the 
specialized  program  makes  it  much  easier  to  produce  mathematical  models  of  sufficient 
realism.  Also,  when  writing  general-purpose  models  of  a  standard  unit  operations,  often  a 
great  deal  of  time  is  spent  preserving  a  degree  of  generality  and  reliability  in  performance 
which  may  not  be  relevant  to  a  particular  problems  [6].  This  type  of  programs  utilizes  two 
different  computational  sequences. 

3.1.1.  Sequential  marching  technique 

In  this  approach,  each  model  equation  is  handled  in  an  order  that  allows  it  to  calculate 
sequentially  the  unknown  variables.  Using  this  technique,  the  mathematical  model 
equations  representing  heat  pump  system  with  refrigerant-filled  collectors  in  space  heating 
applications  was  solved  by  O’Dell  et  al.  [7].  The  system  consists  of  conventional  air 
source,  heat  pumps,  and  active  liquid  solar  collector  units.  Equations  representing  these 
units  are  solved  sequentially.  An  analytical  model  for  a  self-pumping,  boiling  collector 
system  was  developed  and  incorporated  into  a  daily  simulation  by  Walker  and  Davidson 
[8].  The  model  simulates  system  operation  in  three  distinct  isothermal  passes: 
pressurizing,  pumping,  and  running  (heating).  The  conservation  of  mass  and  energy 
equations  for  each  model  are  combined  and  solved  sequentially  and  iteratively.  A  new 
inexpensive  type  of  solar  pond,  mainly  for  space  heating  was  proposed  by  Sokolov  and 
Arbel  [9].  The  system  was  divided  into  four  major  components:  collector,  storage,  ground 
and  climatic  (radiation  and  temperature)  conditions.  These  unit  mode  equations  are  solved 
sequentially  for  predicting  their  system  performance.  El-Nashar  [10]  has  described  the 
mathematical  model  used  to  simulate  the  operational  performance  of  solar  desalination 
plant  which  utilizes  evacuated  tube,  flat-plate  collectors,  and  multi-effect  stack  type 
evaporators,  and  thermally  stratified  heat  storage  tank.  The  major  components  of  the 
process  are  modeled  and  solved  iteratively  in  a  sequence  way.  An  algorithm  based  on  daily 
hours  of  bright  sunshine  was  developed  by  Nguyen  and  Galanis  [11]  The  equations 
representing  this  algorithm  were  solved  sequentially  to  evaluate  the  performance  of  active 
solar  heating  systems  for  different  values  of  system  parameters  (collector  area  and 
efficiency  curves,  orientation  and  inclination,  storage  capacity,  heating  load).  SOLAV 
program  for  numerical  calculation  of  a  solar  heating  system  for  an  individual  building  was 
described  by  Shafeev  [12].  SOLAV  is  one-off  program.  It  consists  of  three  basic  programs 
for  generation  of  climatic  actions  analysis  of  the  solar  heating  system,  with  subprograms 
for  individual  components  of  the  solar  heating  systems  and  determination  of  the  building’s 
heat  demand.  AIRSYS  computer  program  was  also  suggested  by  Shafeev  [13]  for 
designing  active  air  systems  of  solar  heat  supply  for  individual  building.  The  program  was 
used  for  investigation  and  preliminary  estimation  of  the  heat  and  energy  parameters  of 
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such  process.  Equations  representing  a  thermosyphanic  solar  energy  water  heater,  the  heat 
transfer  fluid  circulates  from  the  solar  energy  collector  to  hot  water  store  due  to  the  action 
of  buoyancy  forces,  are  solved  sequentially  by  Norton  and  Probert  [14].  A  methodology 
for  the  constriction  of  a  mathematical  model  of  a  forced  convection  solar  cabinet  dryer  was 
demonstrated  by  Chirarattananon  et  al.  [15].  The  dryer  was  divided  into  three  control 
volumes.  Each  control  volume  model  comprises  only  the  variables  directly  involved  in  the 
energy  and  mass  balance  relationships.  These  mode  equations  are  solved  sequentially  for 
the  state  variables  of  the  control  volumes.  The  concept  of  a  portable  solar  oven  using  a 
vapor-tight  pot  was  introduced  and  theoretically  investigated  by  Ibrahim  and  Khalifa  [16]. 
The  mathematical  model  involved  uses  the  lumped  analysis  approach  where  the  whole 
system  is  broken  up  into  several  elements,  each  of  which  is  treated  as  a  lumped  system  by 
itself.  The  governing  equations  of  these  elements  resulting  from  the  heat  balance  equations 
are  solved  on  the  computer  using  the  fourth-order  Runge-Kutta  method.  A  general  design 
procedure  for  solar-assisted  series  heat  pump  systems  used  for  space  and  process  water 
heating  was  presented  by  Svard  et  al.  [17].  Equations  representing  the  process  units  are 
solved  sequentially.  An  algorithm  to  evaluate  the  performance  of  the  water-cooled  pump 
for  any  location,  given  the  lift  and  collector  area  was  presented  by  Sudhakar  et  al.  [18]. 
The  computer  program  TRNSYS  [4]  was  used  by  Jesch  and  Braun  [19]  to  simulate  the 
performance  of  both  the  fixed  and  variable  volume  thermal  storages.  Two  differential 
equations  describing  the  rate  of  change  of  mass  and  internal  energy  were  given.  These  two 
simultaneous  differential  equations  were  solved  for  each  simulation  time  step.  A  detailed 
computer  simulation  model  of  an  ammonia-water  absorption  cycle  solar  heat  pump  with 
integral  refrigerant  storage  was  reported  by  Kaushik  et  al.  [20].  A  computer  model  based 
on  the  usual  heat  and  mass  balances  for  different  units  in  the  system  was  developed.  The 
unit  models  are  solved  in  sequence  using  one-off  program. 

Elsayed  et  al.  [3]  present  a  textbook  for  design  and  simulation  of  solar  thermal  systems 
using  mathematical  modeling.  Large  number  of  solar  heating  systems  was  considered  in 
this  book.  Mathematical  equations  to  predict  the  radiation  properties  of  multi-layer 
partially  transparent  sheets  with  or  without  an  absorber  plate  are  given.  A  program  to 
calculate  the  radiation  properties  for  several  geometries  of  a  stack  of  similar  or  non-similar 
transparent  sheets  with  or  without  an  absorber  plate  at  any  incidence  angle  was  given.  This 
program  can  also  perform  design  and  optimization  of  flat-plate  collectors  and  solar 
concentrators.  A  computer  program  was  written  by  Elsayed  et  al.  to  design  and  optimize 
the  design  parameters  of  the  flat  collectors,  and  to  calculate  various  factors  affecting  their 
performance.  A  computer  program  to  design  and  optimize  the  design  parameters  of  the 
flat-plate  collectors  was  illustrated.  A  computer  program  was  given  to  assist  the  designer 
in  predicting  the  effect  of  various  designs  and  operating  conditions  on  the  performance  of 
different  solar  desalination  systems.  The  analysis  of  both  the  single  effect  and  the  multi¬ 
effect  were  presented.  Mathematical  simulations  of  various  types  of  sensible  heat  storage 
were  carried  out  for  low-,  intermediate-,  and  high-temperature  applications.  In  particular, 
mathematical  simulations  are  presented  for  the  transient  performance  of  the  mixed  liquid 
storage,  underground  liquid  storage,  and  stratified  liquid  storage.  A  computer  program  was 
also  presented  to  predict  the  performance  of  the  absorption-cooling  machine. 
Mathematical  model  was  given  for  a  Rankin  cycle  selection  of  the  various  design 
parameters.  The  performance  of  the  solar-operated  Rankin  cycle  in  various  operating 
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conditions  was  considered.  Analysis  is  also  presented  to  select  the  optimum  collector 
temperature  for  solar-operated  power  cycles.  A  solar  thermal  desalination  system  with 
heat  recovery  was  recently  modeled  by  Klemens  et  al.  [21].  Numerical  simulation  results 
for  the  process  were  obtained  by  solving  energy  and  mass  balance  equations  for  each  unit 
sequentially.  The  system  units  are  a  solar  collector  and  desalination  tower.  The  modeling 
and  simulation  of  a  solar  water  heating  system  using  time  marching  model  was  presented 
by  Bojic  et  al.  [22].  The  solar  water  heater  consists  of  a  flat-plate  solar  collector,  a  water 
storage  tank,  an  electric  heater,  and  a  water-mixing  device.  A  mathematical  model  of  a 
solar  desalination  unit  with  humidification  and  dehumidification  processes  was  presented 
by  Dai  et  al.  [23]  and  Nafey  et  al.  [24].  Different  mathematical  models  were  proposed  to 
describe  the  simultaneous  heat  and  mass  transfer  process  inside  solar  stills  [25].  The 
process  governing  equations,  including  mass,  momentum,  and  energy  conversations  were 
solved  sequentially  using  one-off  program. 

3.1.2.  Simultaneous  marching  technique 

In  this  approach,  the  considered  processing  unit  or  the  entire  process  is  mathematically 
modeled  in  terms  of  independent  and  dependent  variables.  Then  after  specifications  of  the 
independent  variables,  one  set  of  simultaneous  equations  is  obtained  [6,26,27].  This 
system  of  equations  is  then  solved  by  one  of  the  well-known  numerical  methods  for 
solution  of  systems  of  linear  and  non-linear  equations  [28,29],  to  obtain  the  process 
simulation.  Using  this  technique,  a  space  solar  water  heating  system  was  solved  by 
Dakhoul  et  al.  [30].  The  system  consists  of  four  units,  a  collector,  heat  exchanger,  storage 
tank,  and  the  load  area.  A  set  of  four  non-linear  algebraic  equations  in  four  unknown 
temperatures  was  obtained.  This  set  was  simultaneously  solved  by  Newton’s  iterative 
method.  Adhikari  et  al.  [31]  have  presented  a  computer  simulation  model  for  studying  the 
steady  state  performance  of  a  multi-stage  stacked  tray  solar  still.  The  energy  balance 
equations  for  various  components  of  the  system  were  considered.  A  set  of  simultaneous 
linear  equations  was  obtained  for  water  temperatures  in  various  trays.  This  set  of  equations 
was  cast  in  a  matrix  form  and  solved  iteratively.  A  simplified  simulation  method  for  solar 
thermosyphon  collectors  was  developed  by  Huang  and  Hsieh  [32].  The  process  was 
divided  into  four  parts:  absorber,  storage  tank,  riser,  and  down  comer.  A  set  of  mass  and 
energy  balance  equations  for  these  subsystems  were  simultaneously  solved  by  computer  to 
give  temperature  distribution  and  flow  rate.  The  heat  transfer  analysis  of  an  overlapped 
glass  plate  air  heater  unit  was  presented  by  Selcuk  [33].  Steady  and  unsteady  state 
performances  were  considered.  Eight  simultaneous  differential  equations  were  manipu¬ 
lated  by  finite  difference  technique,  and  in  a  matrix  form.  A  multiple  layer  solar  collector 
was  proposed  by  Kenna  [34].  This  process  was  modeled  by  a  number  of  differential 
equations.  These  equations  were  linearized  and  solved  iteratively.  The  model  was  used  to 
identify  the  most  important  design  parameters  of  this  collector  and  compare  its 
performance  with  flat-plate  collector  characteristics.  A  phase  change  and  sensible  heat 
storage  was  modeled  by  Ghoneim  [35].  Enthalpy  equations  for  a  specified  element,  as  well 
as  the  energy  balance  equation  for  the  fluid  element  were  solved  iteratively  by  Gauss- 
Seidal  method.  A  mathematical  model  derived  from  the  basic  laws  governing  thermal 
energy  exchange  between  surfaces  was  presented  by  Zhao  et  al.  [36]  for  the  transient 
simulation  of  flat-plate  collectors.  A  computer  program  was  coded  in  BASIC  for 
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the  developed  model.  The  concept  of  subdivisions  was  introduced  in  the  development  of 
the  governing  equations  for  an  air  type  solar  collector.  Each  sub  region  set  of  simultaneous 
equations  was  solved  by  Gaussian  elimination  method.  The  subsequent  subdivision 
models  are  solved  iteratively  in  sequence.  Choudhury  and  Garg  [37]  have  developed  a 
mathematical  model  for  jet  plate  air  heater  unit.  A  set  of  energy  balance  equations  for  the 
cover,  absorber,  jet  plate,  back  plate,  air  in  the  passage  between  back  plate  and  air  between 
absorber  and  back  plate  was  solved  iteratively  for  outlet  temperature.  Using  a  program 
specifically  written  for  a  storage  tank  simulation,  the  finite  difference  technique  was  used 
by  El-Nashar  and  Qambiyeh  [38].  A  computer  simulation  of  a  regenerating-type  solar 
collector  using  simultaneous  approach  was  presented  by  Johannsen  [39].  The  collector 
was  used  for  the  regeneration  of  a  liquid  desiccant  and  forms  part  of  an  open-cycle  solar 
air-conditioning  system.  Heat  and  mass  balance  equations  for  the  process  components 
form  a  set  of  13  equations  for  13  unknowns.  A  numerical  model  for  a  single  stage 
humidification-dehumidification  solar  desalination  process  was  presented  by  Fath  and 
Ghazy  [40].  The  desalination  process  consists  of  a  solar  air  heater,  humidifier, 
dehumidifier  and  air-driving  component  (fan).  The  energy  balance  equations  for  the 
process  units  were  coded  by  an  one-off  program  and  solved  simultaneously.  Modeling  and 
performance  analysis  of  a  single  basin  solar  still  with  the  entering  brine  flowing  between  a 
double  glass  glazing  were  investigated  by  Mousa  et  al.  [41].  A  set  of  non-linear  ordinary 
differential  equations  representing  the  process  is  solved  numerically  with  an  adjustable 
time  step  procedure.  An  innovated  water  desalination  system  using  low-grade  solar  heat 
was  studied  recently  by  Al-Kharabsheh  and  Yogi  [42].  The  system  performance  was 
simulated  by  solving  mass,  energy,  and  salt  balance  equations  simultaneously. 

3.2.  General-purpose  programs  (modular  programs) 

Analysis  of  a  number  of  projects  have  indicated  that  the  savings  which  result  from  the 
choice  of  an  optimal  process  configuration  are  far  greater  than  those  resulting  form 
detailed  optimization  of  the  operating  conditions  of  a  plant  of  non-optimal  layout  [43].  The 
choice  of  the  optimal  configuration  will  involve  optimization  of  each  of  the  alternatives. 
This  would  be  accomplished  usually  by  means  of  multiple  case  studies,  which  can  be 
easily  performed  by  a  flexible  general-purpose  program.  In  other  word,  in  order  that  the 
computer  program  can  be  applied  widely  for  different  solar  heating  processes,  with 
different  configurations  and  different  operating  conditions  it  must  be  flexible  and  general. 
To  meet  these  requirements,  the  modular  structure  should  be  adopted  [44]. 

In  this  structure  the  mathematical  model  is  usually  formulated  in  terms  of  sets  of 
equations  relating  to  the  units  of  the  process.  Each  of  these  sets  of  equations  may  be 
regarded  as  an  independent  and  self-standing  module.  This  is  the  basis  of  the  modular  flow 
sheeting  approach.  Once  a  module  has  been  developed,  it  can  be  used  in  any  number  of 
processes  requiring  such  a  module.  Modules  represent  a  set  of  ‘building  blocks’  which  are 
selected  and  assembled  by  the  user  to  produce  a  model  of  the  complete  process  he  wishes 
to  study.  The  flow  sheeting  approach  had  been  used  in  the  design  of  chemical  processes  for 
a  long  time.  Similar  approach  was  in  fact  used  in  other  fields,  such  as  electrical 
engineering,  thermal  design  systems,  but  these  were  referred  to  by  other  names.  Excellent 
reviews  of  the  early  approaches  to  flow  sheeting  were  given  by  Westerberg  et  al.  [45] 
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and  Wells  and  Robson  [46].  Flow  sheeting  programs  have  manipulated  two  computational 
sequences. 

3.2.1.  Sequential  modular  approach 

In  the  sequential  modular  approach,  the  model  equations  are  handled  using  a  library  of 
‘modules’  or  subroutines.  These  modules  are  able  to  compute  output  streams  for  given 
input  streams  for  each  equipment  in  the  simulated  process.  There  is  then  an  executive, 
which  uses  the  process  flow  structure  to  link  the  subroutines  in  order  to  perform  the 
required  computation  for  the  whole  process.  Specification  such  as  input  stream  data  and 
equipment  parameters  are  easily  handled  by  passing  the  specified  values  directly  to  the 
proper  unit  modules.  However,  other  specifications,  on  output  streams  for  instance,  may 
not  accept  the  specifications.  Other  specifications  that  cannot  be  input  directly  to  a  module 
are  referred  to  as  design  specifications. 

SYSOPT  modular  program  was  described  by  Young  [47].  This  program  uses  linear 
programming  technique  to  optimize  the  subsystem  sizing,  performance,  and  economics  of 
integrated  energy  systems.  The  program  modular  stricture  permits  the  designer  to  readily 
interconnect  solar  subsystem  models  into  alternative  configurations  of  combined  solar  and 
conventional  technologies.  Based  on  these  interconnections  and  specifications,  the 
program  formulates  hourly  energy  /mass  flow  balances  at  interconnection  points  and 
optimizes  both  sizing  and  time-dependent  dispatch  of  the  various  subsystems.  The  calling 
of  the  available  subsystem  modules  in  the  optimization  program  SYSOPT  includes  solar 
flat-plate  collectors,  photovoltaic  collectors,  wind  energy  systems,  electric  storage, 
thermal  storage,  turbine  generators,  power  conditioners  and  inverters,  resistance  heating, 
thermal  and  electric  load  modules,  fossil  boiler,  diesel  generators,  utility  grid  backup  and 
feedback  and  various  process  modules  such  as  water  desalination.  Additional  subsystems 
models  can  be  readily  formulated  within  the  overall  modular  structure.  Each  subsystem 
modules  is  specified  by  a  single  subroutine,  which  returns  the  appropriate  coefficients  for 
each  subsystem  variable. 

TRNSYS  (Transient  Simulation  Program)  is  a  thermal  process  simulation  program.  A 
list  of  different  components  in  the  program  library  was  presented  by  Duffle  and  Beckman 
[4].  The  unit  mathematical  model  was  represented  by  a  number  of  algebraic  and/or 
differential  equations.  Three  different  integration  techniques  were  used  by  TRNSYS 
program.  A  number  of  simulation  studies  for  different  solar  heating  processes  were 
reported  by  Duffle.  Duffle  hinted  five  modular  computer  programs.  These  are:  WATSUN 
program,  which  was  developed  at  the  university  of  Waterloo  in  Canada.  This  program  can 
manipulate  a  large  number  of  unit  modules.  MINSUN  program  is  specifically  designed  for 
simulation  of  systems  with  seasonal  storage  where  daily  temperature  changes  are  small; 
one-day  time  step  are  used.  G3  is  a  simplified  program.  It  is  based  on  daily  input-output 
diagrams.  This  program  is  designed  for  fast  and  simple  evaluation  of  active  systems.  ISFH 
is  a  versatile  program.  It  includes  a  means  of  generating  a  ‘synthetic  climate’  from 
monthly  average  weather  data.  EMPG2  (European  Modeling  Group  Program  2)  is  a 
modular  program.  It  has  the  capability  of  handling  a  wide  verity  of  transient  thermal 
systems,  system  control  strategies,  and  load  types.  Also,  several  programs  for  different 
solar  heating  processes  have  been  developed.  These  programs  were  used  as  modules  for 
TRNSYS  modular  program.  For  example,  a  simulation  model  for  fresh  water  production 


A.S.  Nafey  /  Renewable  and  Sustainable  Energy  Reviews  9  (2005)  576-591 


585 


from  a  desalination  system  integrated  in  a  greenhouse  roof  was  developed  by  Chaibi  [48]. 
This  module  was  based  on  a  set  of  heat  balance  equations  including  material  layers  and  a 
brackish  water  layers.  The  equations  were  solved  with  the  more  general  PC  program, 
Engineering  Equation  Solver.  For  calculation  of  solar  irradiation,  the  simulation  code 
TRNSYS  was  used.  A  preliminary  study  of  a  solar-heated  low  temperature  space  heating 
system  with  seasonal  storage  in  the  ground  was  performed  [49].  The  system  performance 
was  evaluated  using  simulation  module  (DST)  [50]  for  duct  ground  heat  storage  model. 
DST  program  was  used  as  a  module  for  the  modular  program  TRNSYS. 

The  following  are  the  advantages  of  sequential  modular  approach  [51]:  on  the  module 
level,  one  or  more  specialized  algorithms  can  be  used  with  each  module  to  solve  the  model 
equations  thus  the  module  calculations  can  be  very  efficient  and  robust.  The  engineer 
easily  understands  the  sequential  modular  programs,  because  the  information  flow  in  the 
program  closely  resembles  the  material  flow  in  the  process.  Also,  the  information  flow 
structure  makes  error  checking  fairly  easy. 

The  disadvantages  are  [6] 

•  The  sophisticated  handling  of  design  specifications. 

•  The  presence  of  multiple  nested  iteration  loops. 

•  Not  well  suited  to  process  optimization. 

The  handling  of  design  specifications,  typically  equations  of  the  form  x  =  constant,  by 
introducing  additional  iteration  loops  is  obviously  an  inefficient  way  to  handle  such  simple 
equations.  The  second  problem  can  be  explained  by  noting  that  the  accelerated  subroutine 
loops  used  to  converge  the  tear  stream,  together  with  any  control  loops  used  to  handle 
design  specifications,  may  represent  the  outmost  loops  in  a  hierarchy  of  nested  iteration 
loops.  Within  these  loops  are  many  iteration  loops  that  may  be  needed  within  the  process 
unit  modules,  and  at  the  innermost  level  any  loops  that  may  exist  within  the  physical 
property  subroutines  called  by  the  unit  modules.  Problems  involving  large  systems  of 
nested  iteration  loops  would  be  very  expensive  and  time-consuming  to  solve.  Finally,  the 
optimization  problem  cannot  be  handled  easily  by  the  sequential  modular  approach, 
because  adding  more  outer  iteration  loop  for  optimization  will  complicate  the  numerical 
solution. 

Because  of  the  fundamental  problems  described  above,  there  has  been  recent  interest  in 
an  alternative  approach  to  process  simulation  and  optimization.  This  approach  namely  the 
simultaneous  modular  approach  is  discussed  below. 

3.2.2.  Simultaneous  modular  approach 

In  the  case  of  the  simultaneous  modular  approach,  the  simulation  problem  is  seen  as 
being  one  of  solving  a  large  set  of  linear  and  non-linear  equations.  The  modules  present  the 
executive  with  the  equations  describing  the  unit  process.  The  executive  assembles  all  these 
equations  and  chooses  the  optimum  method  of  solution  for  the  whole  system.  The 
QUASILIN  equation-oriented  flow  sheeting  program  was  described  by  Hutchison  et  al. 
[52].  Its  performance  on  several  simulation,  design,  and  optimization  problems  was 
reported.  Newton’s  method  was  used  for  equation  solving  and  the  Powell  successive 
quadratic  programming  algorithm  for  optimization. 
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Though  the  simultaneous  approach  is  commonly  used  in  chemical  and  thermal 
processes  simulation,  not  many  studies  for  solar  heating  process  simulation  using  this 
technique  have  been  published  so  far. 

On  one  hand  this  approach  provides  flexibility  in  handling  design  constraints  and  it  has 
advantages  in  problems  with  recycle  structure  [1,45,53].  Perhaps  the  most  important 
benefits  of  this  approach  arise  when  performing  optimization  performance  calculations  or 
dynamic  simulation.  Here  the  ability  to  solve  repeatedly  and  quickly  a  large  set  of  non¬ 
linear  algebraic  equations  are  useful. 

On  the  other  hand,  the  disadvantages  are  [54]  great  effort  required  in  the  development 
of  an  intelligent  program  for  coordination.  A  universal  numerical  procedure,  which  cannot 
be  adapted  to  the  structural  properties  of  the  individual  blocks,  is  needed. 


3.3.  Structure  of  a  simulator 

From  the  early  stage  of  having  special  purpose  (one-off)  program  for  simulating  one 
process  or  one  unit,  the  art  of  solar  heating  process  simulation  has  developed  into  an 
advanced  form,  where  proficient  programs  can  simulate  a  large  number  of  processes  and 
configurations. 

The  main  components  of  a  simulation  package  are  shown  in  Fig.  3.  The  package  can  be 
a  fixed  structure  or  a  flexible  structure.  Programs  with  rigid  structure  are  simple  and 
straightforward.  All  the  user  has  to  supply,  is  the  data,  and  the  executive  handles  the 
program  in  the  same  way,  irrespective  of  the  sort  of  the  process  simulated.  Unfortunately, 
all  unit  modules  have  to  be  loaded  into  the  computer  main  memory  whether  or  not  many  of 
the  routines  are  relevant  for  the  process  being  simulated.  As  these  programs  grow  to 
include  more  unit  modules,  physical,  and  thermodynamic  property  calculations  and  wider 
output  options,  the  complexity  increases. 


Fig.  3.  Structure  of  a  modular  program. 
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Fortunately,  most  modern  simulators  use  a  flexible  or  variable  structure.  In  the 
flexible  structure,  a  problem-oriented  language  (POL)  is  used.  The  process  topology 
can  be  explained  by  sentences.  Only  those  modules  that  are  actually  needed  are 
loaded.  Data  storage  space  may  be  allocated  only  as  required  by  the  size  of  the 
problem.  The  variable  structure  simulator  provides  greater  flexibility  and  modularity. 
From  Fig.  3,  it  is  seen  that  the  structure  of  a  simulator  essentially  consists  of  five 
phases: 

•  input  and  input  verification 

•  calculation 

•  processing 

•  output 

The  input  should  consist  of  process  topology,  information  on  all  input  streams,  design 
parameters.  Process  topology  is  defined  through  POL.  Feed  stream  information  consists  of 
flow  rate,  composition,  temperature,  and  pressure.  The  design  parameters  of  the  units  may 
be  solar  collector  area,  heat  exchanger  area,  etc.  Convergence  criteria  would  normally  be 
the  tolerance  on  the  convergence  of  recycles  calculations.  The  maximum  number  of 
iterations  allowed  may  also  be  specified. 

The  processing  stage  consists  of  generating  the  main  program  in  a  procedure-oriented 
language,  compiling  the  program,  loading  all  the  required  modules  for  the  simulation, 
linking  the  main  program  and  all  subprograms,  and  allocated  storage  prior  to  the  execution. 
Most  simulators  have  extensive  input  data  checking  facilities  as  part  of  the  processing 
phase. 

The  calculation  stage  requires  unit  module  algorithms  or  equation  sets,  physical 
properties  data,  thermodynamic  property  estimation  routines,  numerical  methods  for 
convergence  promotion  and  any  user-supplied  algorithms.  In  almost  all  simulators,  a 
stream  table  is  a  common  feature  of  output.  The  table  consists  of  stream  flow  rate, 
composition,  temperature,  pressure,  and  enthalpy  for  all  streams  in  the  process.  The  post¬ 
processing  phase  would  generally  constitute  cost  estimation  and  economic  analysis. 


4.  Numerical  techniques 

Computer  solutions  of  individual  solar  heating  process  units  are  often  highly  iterative 
[3].  Often  the  models  of  these  units  involving  rather  sophisticated  numerical  methods  for 
rapid  convergence  [55].  By  connecting  these  individual  units  to  form  entire  process,  one  or 
more  recycle  streams  may  be  involved.  When  this  situation  exists,  it  often  becomes  very 
difficult  to  obtain  a  numerical  solution.  If  more  than  one  recycle  loop  is  present,  the 
problem  tends  to  grow  exponentially.  It  is  also  possible  that  the  particular  process  is 
computationally  unstable,  meaning  that  a  solution  cannot  be  obtained  by  conventional 
techniques.  Cavett  [56]  has  analyzed  the  problems  of  simulated  processes  involving 
recycle  loops  in  terms  of  existing  numerical  mathematics  and  attempts  to  apply  more 
sophisticated  techniques  to  their  solution.  Several  methods  are  considered,  including 
successive  substitution,  relaxation,  and  Newton’s  method. 
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Generally,  solar  heating  process  simulation  gives  rise  to  the  need  for  solving  sets  of 
linear  and  non-linear  algebraic  equations.  For  steady-state  calculations,  the  equations  are 
naturally  algebraic,  while  dynamic  problems  involving  mixed  systems  have  algebraic  and 
differential  equations  whose  solution  in  turn  involves  iterative  solution  of  algebraic 
systems.  This  problem  is  most  effectively  dealt  with  by  either  Newton’s  Raphson  (NR) 
method  or  a  quasi-Newton  method.  However,  NR  method  is  usually  a  slower  method  than 
a  specially  tailored  numerical  scheme  for  solving  the  equations  describing  a  particular 
unit.  A  wide  variety  of  iterative  solution  procedures  for  solving  linear  and  non-linear 
equations  are  adopted  for  solving  solar  heating  processes  [3]. 


5.  Program  evaluation  (numerical  results  validity) 

The  question  considered  in  this  section  is  simple,  but  fundamental.  How  do  we  know 
that  the  results  obtained  numerically  are  correct? 

To  build  confidence  in  large-scale  computer  programs  for  solar  heating  systems 
simulation,  the  programs  must  be  reliable,  free  of  bugs,  and  a  large  number  of  different 
problems  must  be  solved  by  it.  In  assessing  the  reliability  of  the  final  results,  the  accuracy 
and  precision  should  be  considered  [43].  Accuracy  may  be  expressed  as  how  closely  the 
experimental  results  agree  with  the  true  or  most  probable  values.  In  contrast,  precision 
may  be  defined  as  how  closely  a  set  of  measurements  of  the  same  quantity  agrees  with 
each  other.  Thus,  accuracy  expresses  the  correctness  of  the  measurement  whereas 
precision  describes  the  predictability  of  the  measurements. 

The  validation  activity  basically  differs  from  classical  comparisons  between 
experimental  data  and  theoretically  predicted  values,  because  it  accounts  not  only  for 
the  measuring  errors  to  be  assigned  to  the  experimental  data,  but  also  for  the  uncertainty  to 
be  associated  with  the  calculated  values  as  a  result  of  errors  on  measured  input  data  and  the 
uncertainty  about  the  correct  value  of  the  various  parameters  describing  the  process. 
Experimental  evaluation  for  some  of  the  published  mathematical  models  for  water 
desalination  by  solar  still  was  presented  by  Abd  El  Kader  et  al.  [25].  The  validity  and 
applicability  of  these  models  under  real  operating  conditions  were  examined. 


6.  Future  development 

The  availability  of  inexpensive,  yet  powerful  personal  computers  has  now  opened  new 
applications  for  the  development  of  programs  adapted  to  these  personal  computers.  One  of 
the  essential  aspects  for  the  development  of  software  for  personal  computers  is  user 
friendliness.  In  other  words,  the  user  must  be  in  a  position  to  enter  the  process  topology 
and  data  interactively.  Also,  the  user  must  be  able  to  modify/change  any  solar  heating 
process  unit  with  minimum  effort  and  maximum  efficiency  and  flexibility.  This  can  be 
achieved  by  using  graphical  user  interface.  This  approach  is  a  strong  move  from 
the  representation  of  a  system  and  its  performance  as  simple  lists  of  numbers  and  letters, 
the  characteristics  of  the  most  present  programs,  to  graphical  representations  that  can  be 
achieved  by  graphical  interface. 
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The  use  of  exergy  analysis  for  rating  and  optimizing  the  performance  of  energy 
conversion  systems  is  accepted  on  a  worldwide  basis  as  one  of  the  most  reliable  tool 
available  for  this  purpose  [5,57].  Using  this  analysis,  many  solar  heating  system  units  have 
been  developed  [59-61].  However,  the  exergy  concept  has  not  yet  been  used  in  a  general 
modular  program.  By  considering  the  entropy  generation  as  well  as  the  material  and 
energy  balances  equations  a  realistic  model  can  be  obtained.  Each  individual  component 
or  assembly  can  then  be  optimized,  by  minimizing  the  respective  entropy  generation. 


7.  Summary 

It  is  evident  from  the  foregoing  discussion  that  simulation  is  a  very  important  tool  for 
performance,  design,  and  optimization  of  solar  heating  systems.  These  systems  which  are 
governed  by  the  principles  of  heat  transfer,  thermodynamics,  fluid  mechanics,  and  mass 
transfer,  arise  in  a  wide  range  of  solar  engineering  applications.  The  mathematical  models 
representing  these  applications  may  comprise  a  large  system  of  equations  caramelized  by 
complexity,  non-linearity,  and  implicitly.  Many  computer  programs  for  different  solar 
heating  processes  are  reported.  These  programs  are  classified  in  this  work  into  two  types. 
First,  special  purpose  programs  (one-off  programs)  and  second,  general-purpose  programs 
(or  modular  programs).  Various  aspects  of  these  programs  such  as  the  structure, 
computational  technique,  advantages,  and  disadvantages  of  these  programs  are  illustrated. 
The  basic  considerations  involved  in  numerical  modeling,  particularly  those  concerned 
with  accuracy  and  validation,  are  discussed. 

The  thermal  engineer  involved  with  process  simulation  is  working  in  a  very  dynamic 
environment.  Two  points  for  future  development  are  pointed.  They  are  using  the  exergy 
concept  in  modeling  a  general  modular  program,  and  using  the  new  feature  of  the  visual 
programming  in  building  flexible,  and  easy  to  manage,  programs  with  comfortable 
interface. 
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